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Abstract 

In the context of Landau theory and its field theoretical refine- 
ments, interfaces between coexisting phases are described by intrinsic 
profiles. These intrinsic interface profiles, however, are neither di- 
rectly accessible by experiment nor by computer simulation as they 
are broadened by long-wavelength capillary waves. In this paper 
we study the separation of the small scale intrinsic structure from 
the large scale capillary wave fiuctuations in the Monte Carlo sim- 
ulated three-dimensional Ising model. To this purpose, a blocking 
procedure is applied, using the block size as a variable cutoff, and a 
translationally invariant method to determine the interface position 
of strongly fluctuating profiles on small length scales is introduced. 
While the capillary wave picture is confirmed on large length scales 
and its limit of validity is estimated, an intrinsic regime is, contrary 
to expectations, not observed. 
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1 Introduction 

Under suitable conditions systems of statistical mechanics can exhibit co- 
existence of different phases which are separated spatially by interfaces. 
Structure and properties of such interfaces have been investigated both 
experimentally and theoretically, and still offer problems for current re- 
search. Near critical points several interfacial properties are expected to 
show universal behaviour. 

Since the 19th century, there has been continuous interest in the the- 
oretical description of interfaces ^ El El E]- Traditionally, the interface 
has been described by mean field theories and their refinements El El 
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as an "intrinsic" continuous profile with a width proportional to the bulk 
correlation length. 

The existence of the interface, however, breaks the translational invari- 
ance of the system, resulting in long-wavelength capillary wave fluctua- 
tions of the interface position as Goldstone- modes. These fluctuations are 
neglected by mean field theory which assumes a flat interface, but they 
strongly influence the interfacial properties. According to capillary wave 
theory jHUD, the interface width, for example, depends logarithmically on 
the system size and diverges in the limit of an infinite system. 

Mean field and capillary wave theory can be viewed as a small-scale and 
large-scale description of the interface, respectively, and then be combined 
in the "convolution approximation" 0m . In this picture the interface is 
described as an intrinsic profile which is centered around a two-dimensional 
surface subject to capillary wave fluctuations. 

While capillary wave behaviour has been observed in many systems (e.g. 
experimentally on liquid- vapor interfaces jJOl and polymer mixtures ^I] , in 
simulations of polymer mixtures and of the Ising model \V3l 1141 [TE\ ). 
the mean field behaviour is difficult to access as the intrinsic profile - if 
it has a well-defined meaning at all - is broadened by the capillary wave 
fluctuations. 

The concept of intrinsic versus large-scale interface structure has so 
far not been defined unambiguously outside a given theory. There is no 
clear "snapshot" definition of an interface, which always shows overhangs, 
islands, clusters which may be as well assigned to the phase-separating 
interface region as to the fluctuation structure of the neighbouring homo- 
geneous phases. 

In this paper the attempt is made to separate the intrinsic structure 
from the influence of the capillary waves via a blocking procedure reminis- 
cent of the Kadanoff block spin method of renormalization group theory. 
The system of size L x L x D with an interface perpendicular to the in- 
direction is divided into blocks of size B x B x D. The different length 
scales are separated by using the block size B as a variable cutoff. To cal- 
culate the interface position of the strongly fluctuating interface profiles on 
small length scales, we implement a method which respects translational 
invariance (in contrast to various other approaches in the literature). In 
this way, the local interface position, the interface profile and the interface 
width can be studied on the different length scales, thereby allowing a di- 
rect test of the convolution approximation and of the limits of validity of 
the capillary wave theory and of the mean field theory. It is found that on 
large length scales, the capillary wave picture is consistent with the data. 

On small length scales, however, contrary to expectations we cannot 
identify an intrinsic interface profile with a width of the order of the corre- 
lation length. For the case of the Ising model under study the concept of 
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the intrinsic profile is thus called into question. 

An approach based on blocking has already been used by Weeks [TB] 
for liquid-vapour interfaces in order to derive the capillary wave model 
from a microscopic model Hamiltonian, thereby reconciling mean field and 
capillary wave theory as small-scale and large-scale interface descriptions, 
respectively. In the context of Monte Carlo simulations a procedure simi- 
lar to the one used in this paper has been applied by Werner et al. ^21 to 
homopolymer interfaces, finding consistency of their data with the convo- 
lution approximation. 

2 Theoretical description of the interface 
profile 

2.1 Intrinsic interface profile 

Interfaces can be described in mean field theory by a continuous profile 4>{x) 
representing the difference between the concentrations of the two coexist- 
ing phases. In the framework of Landau theory and its field theoretical 
extensions the profile function plays the role of a local order parameter. 
Near the critical temperature the free energy density is written as 17,: 

c^^idcj^ixyr + viHx)) (1) 

with the double-well potential 

¥{<!>) = A(02 - 1)2 . (2) 

The minima of this potential correspond to the two homogeneous phases, 
whose densities have been normalized to ±1 for simplicity. 

Minimization of the free energy density with boundary conditions 
appropriate for an interface perpendicular to the z-axis leads to the typical 
hyperbolic tangent mean field profile 

=tanh(^^(z-/i)^ (3) 

with a width proportional to the mean field correlation length which 
diverges near the critical point with exponent 1/2. The parameter h spec- 
ifies the location of the interface which is arbitrary due to translational 
invariance. 

Corrections due to fluctuations of the order parameter field can be cal- 
culated systematically in renormalized perturbation theory. In the "lo- 
cal potential approximation" only corrections to the local potential V{(j)) 
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are taken into account, whereas contributions involving higher powers of 
derivatives are neglected. The resulting interface profile is of the form 
(p{z) = f{{z — /i)/2^), where the width is now proportional to the physical 
correlation length ^, which diverges near the critical point with exponent 
ly. This corresponds to the scaling form proposed by Fisk and Widom [Jj. 
To lowest order the profile function / is equal to the mean field profile 
function, / = tanh. In higher orders, corrections to the potential V{(j)) due 
to perturbative loop contributions modify the mean field profile |18| . We 
do, however, not consider them as they are numerically small and can be 
neglected in our context. 
The interface profile 



in lowest order of the local potential approximation is to be distinguished 
from the mean field profile. It represents a refinement of Landau theory 
that is consistent with scaling [T] and holds near the critical point. 

In the local potential approximation long wavelength fluctuations are 
not fully taken into account. Consequently the profile does not contain 
the effects of capillary waves, which are discussed in the following section. 
The function above thus represents an intrinsic interface profile. 

2.2 Capillary wave theory 

Translational invariance is broken by the presence of the interface. The 
Goldstone bosons associated with this broken symmetry are long wave- 
length excitations of the interface position, which have vanishing energy 
cost in the infinite wavelength limit 19 . These capillary waves strongly 
influence interfacial properties like the interface width but are neglected by 
mean field theory which assumes a flat interface. In capillary wave theory, 
as introduced by Buff et al. bubbles, overhangs and any continuous den- 
sity variations of the interface are neglected. The local interface position 
is parameterized by a single- valued function h(x,y), describing the inter- 
face as a fluctuating membrane. The free energy cost of capillary waves is 
basically due to the increase in interface area against the reduced interface 
tension a. Thus the capillary wave free energy (multiplied by the inverse 
temperature (3) is [2] 



where |V/i| ^ 1, i.e. long wavelength and small amplitude, has been as- 
sumed. The Gaussian nature of the capillary wave model allows analytic 





(5) 
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calculation of thermal averages. The local interface position is distributed 
according to a Gaussian distribution 



P{h) = {S{h{x, y) - h)) = 0-'^'/'^' (6) 

V27r s 



with variance 



s' = {h^) = ^ ldq^- = ^\r.^. (7) 
iTTcr J q lira Be 



■'cw 



In order to avoid divergence of the integral, a lower cutoff 2'k/L and an 
upper cutoff 27r/i?cw have been introduced. Here L is the system size, 
which determines the maximal allowed wavelength and provides a natural 
lower cutoff for the capillary wave spectrum. The cutoff length i3cw, which 
is so far arbitrary, sets the scale below which capillary wave theory ceases 
to be valid. As there should be no capillary waves with wavelength smaller 
than the intrinsic width of the interface, Sew should be of the order of 
the correlation length. The factor 27r in the lower cutoff is due to periodic 
boundary conditions and has also been introduced in the upper cutoff for 
convenience. The phenomenological capillary wave model thus requires two 
inputs: the macroscopic interface tension a and the cutoff i?cw While the 
interface tension near the critical temperature is known from scaling 0], 
the cutoff J5cw is unknown. 

In the capillary wave picture the instantaneous interface profile is a 
sharp step function between the two phases. Nevertheless, in the thermal 
average the capillary wave fluctuations produce a continuous density profile 

p{z) = Jdh sgn{z - h) P{h) = cii (^-^^ (8) 

with a finite width whose square is proportional to the variance and 
thus to the logarithm of the system size L. In other words, the apparent 
width of the interface depends on the length scale on which the interface 
is studied and diverges in the thermodynamic limit. 



2.3 Convolution approximation 

The two different descriptions of the interface profile - intrinsic profile and 
capillary wave theory - can be reconciled by adopting the view that the 
intrinsic profile describes the interface on a microscopic scale of the order 
of the correlation length, while capillary wave theory describes the macro- 
scopic interface fluctuations of wavelengths much larger than the correlation 
length. Neglecting the coupling of the capillary wave fluctuations to the 
intrinsic interfacial structure, the interface can be viewed as a fluctuating 
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capillary wave theory membrane h(x, y), to which an intrinsic profile (f>{z) is 
"attached" . In this way an instantaneous profile (j){z — h{x,y)) is obtained. 
Performing the thermal average over the capillary waves, the "apparent" 
profile is given by the convolution 

m{z)= Jdh(f>{z-h)P{h). (9) 

In this convolution approximation, the intrinsic profile is thus broadened 
by capillary wave fluctuations. The broadening increases with the system 
size according to Eqs. © and {T)) corresponding to the fact that in a larger 
system more capillary waves are allowed. 



2.4 Interface width 

The interface width cannot be defined in a unique way, and various def- 
initions have been used in the literature (see e.g. refs. 0, [S] and ^^)- 
For numerical purposes, a suitable definition of the squared interface width 
is the second moment of a weight function which is proportional to the 
gradient of the interface profile m: 

. I^!±^. (10) 

Jdz m'[z) 

Due to the linearity of this definition and of the convolution approximation 
© , the squared width w"^ of the convolution profile m is simply the sum 
of the intrinsic width squared wf^^. of the intrinsic profile (p, Eq. Q), and 
of the capillary width squared w'^^ = of the capillary wave profile p, 
Eq. ®: 

TT^ 1 L 

This expression shows the broadening of the intrinsic profiles by capillary 
waves. It can be cast into a temperature-independent form by expressing 
all lengths in units of twice the correlation length: 

wL^yY^j =U;int+Wew = Y^+^cwl0g2-g— , (12) 

where L = L/2£^, B^w = Here the slope 

In 2 



^'cw - - 0.265(2) (13) 

contains the the universal constant a^^ ~ 0.1040(8) |20[ I21| . We also con- 
sidered other definitions of the interface width, e.g. with a weight function 
proportional to (to'(z))^, but the results to be discussed below are not 
altered significantly. 
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3 Monte Carlo calculation 



3.1 Monte Carlo setup 

In experiments or Monte Carlo simulations of systems with interfaces the 
observed interface profile and width are the total ones, including the effects 
of the intrinsic structure as well as of the capillary waves. In order to inves- 
tigate the question whether the intrinsic structure can be separated from 
the effects of capillary waves, suitable observables have to be constructed. 
We have studied this question by means of Monte Carlo simulations of 
the three-dimensional Ising model, which is in the same universality class 
as three-dimensional binary mixtures. The results allow to test the pre- 
dictions from capillary wave theory and the validity of the convolution 
approximation. 

In the simulations a cubic lattice of size L x L x D is used, with spins 
s{x, y, z) £ {±1} on each lattice site (x, y, z). The presence of an interface 
is enforced by antiperiodic boundary-conditions in the z-direction, while 
the boundary conditions in the x, y-directions are periodic. The lattice 
sizes are varied from L = 32 to 512 and from _D = 50 to 180. To reduce 
critical slowing down, the Single-Cluster- Wolff algorithm ^] is used to 
allow simulations at reduced temperatures t = {T — Tc)/Tc = —0.05, —0.01 
and —0.004, corresponding to correlation lengths ^ — 1.65, 4.56 and 8.13 
in lattice units, respectively. 

In order to analyze the interface on different length scales a blocking 
procedure is applied. The system is split into columns of block size B x B 
and of length D. The block size B will later be varied to allow for systematic 
analysis. In each block i a local block profile is calculated by averaging over 
the lateral coordinates in the block: 

™*(^) ^ X! s{x,y,z). (14) 

2:,y(iBlock i 

3.2 Local interface position 

Due to strong thermal fluctuations near the critical point the interface 
position for this profile is not well-defined. Various methods have been 
proposed to determine the interface position (see e.g. refs. 22|> [221 and 
|24j ) and applied to total interface profiles, i.e. for the profile l)14|l at B = L. 
These methods, however, are not appropriate for our purpose, because they 
either neglect bulk fiuctuations, whose contributions we want to take into 
account fully, or they suffer from a lack of translational invariance: If the 
interface is, due to interface wandering, not in the middle of the system, 
the fluctuations above and below the interface do not average out, leading 
to a biased value for the interface position. This bias can be neglected 
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for profiles on sufficiently large length scales, i.e. for block profiles with 
sufficiently large B, but not for the strongly fluctuating profiles on smaller 
length scales B. 

To solve this problem, a translationally invariant method to determine 
the interface position is proposed |2S1-^ Because of the periodic/antipe- 
riodic boundary conditions, the Ising lattice is wrapped to a torus. The 
antiperiodicity is implemented by taking the couplings between nearest 
neighbours s(x, y, z) and s{x, y, z + f ) to be negative for a particular value 
oi z = Za- Superficially, the choice of Za appears to break translation in- 
variance along the ^-direction. This is, however, not the case. Under a shift 
of Za the Hamiltonian remains invariant, if at the same time those spins, 
which cross the plane of antiperiodicity, are flipped. Consequently trans- 
lation invariance remains a symmetry. Observables should be constructed 
in such a way that they also respect translation invariance. Translational 
invariance of our method to determine the interface position is achieved 
in the following way. The antiperiodicity point Za is shifted through the 
system successively. The values of the block profiles mi(z), which cross the 
antiperiodicity plane z = Za, change signs appropriately. For a monotone 
interface profile function m(z), all values would get the same sign, if Za 
coincides with the center of the interface. Correspondingly, we define the 
location z of the interface to be that value of Za, which makes the number 
of values of mi{z) having the same sign maximal, i.e. being mostly positive 
or mostly negative. Mathematically this point can be determined as the 
maximum of | '^iC'^)!- ^^^^ way, one assures that the fluctuations 
near the interface determine the precise location of the interface position, 
ft would also be possible to use the minimum of | mi{z)\, which is at- 
tained when the interface is located opposite of Za- fn this case, however, 
fluctuations deep in the bulk phase opposite (on the torus) to the interface 
position would account for the precise location of the interface position, 
which is physically less meaningful. 

This boundary shift method very successfully determines the interface 
positions for strongly fluctuating profiles even on rather small length scales 
B, but of course fails for extremely small length scales B — f,2,3, where 
the profiles mi(z) are still reminiscent of their Ising nature with only two 
allowed values ±1 and thus do not possess a well-defined interface position 
at all. This will result in Ising lattice artefacts for extremely small block 
sizes in the later discussion. But as neither of the interface models discussed 
in section El applies on a microscopic scale, this deficiency of the method 
does not prevent a test of the models. 

^After completion of this work ref. I'ltil has been pointed out to us, where the same 
method is being used. 
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3.3 Interface profile 

In order to obtain the interface profile without fluctuations on length scales 
larger than B, the block profiles mi{z) can now be shifted such that the 
interface position is in the middle of the system. Note that values that pass 
the antiperiodic boundary conditions change sign (spin flip) . To assure that 
the profile values above the interface are predominantly positive, an overall 
sign change is applied if necessary. The shift procedure corresponds to a 
measurement in a moving frame that is wandering with the local interface, 
i.e. to an elimination of the zero mode within the block. One can now 
take the block and Monte Carlo average over the block profiles mi{z) to 
obtain the local interface profile m(z) on the scale of the block length 
B. This profile corresponds to the apparent profile ^ of the convolution 
approximation measured on a length scale of the block size B. 

The interface width of the profile is determined via a discretized gradient 
analysis corresponding to Eq. H1Q|I . Let the z-coordinates take the values 
1,2, . . . , D. An auxiliary coordinate t — z — D/2 is introduced which labels 
the positions between adjacent lattice sites and is centered about the middle 
of the system, assuming the values —D/2 + 1, —D/2 + 1, . . . , D/2 — 1. A 
normalized gradient is then defined by 



with the normalization 



m[t+ — + l\ -m\t+ — 



toU+ — + l -m t+ — 



(15) 



(16) 



2 J V 2 
The interface width is the second moment of this weight function: 



So far each block has been considered as an Ising system of size BxBxD 
in its own right, for which the same measurements as for a total Ising system 
have been made. In this way all fluctuations of wavelengths larger than 
B have been cut off. The opposite approach is to flrst average over each 
block and then consider the resulting Ising system of size L/B x L/B x D, 
thereby "averaging out" all fluctuations of wavelength smaller than B. 
This coarse graining view is adopted in order to obtain a capillary wave 
like description of the interface. In each block i the interface position hi 
is determined by the above described boundary shift procedure and then 
measured relative to the total interface position, i.e. the position of the 
profile rrii for B — L, determined by the same method. The resulting 



9 



height variables for each Monte Carlo configuration form the snapshot of 
a membrane corresponding to the function h(x, y) of the capillary wave 
model, but coarse grained on the length scale B. According to the capillary 
wave model, the distribution P{h) of these height variables hi over all blocks 
and Monte Carlo configurations should be Gaussian, Eq. 

4 Results 

4.1 Confirmation of the capillary wave model 

The measured distributions P{h) of the interface heights a,t t — —0.01 
and L = 128 for various coarse graining lengths B are shown in Fig. ^ 
They display Gaussian-like peaks which become narrower with increasing 
block size, thus qualitatively following the capillary wave model predictions 
Eqs. ® and 0. For B = L the distribution has degenerated to a delta 
peak. The additionally displayed Gaussian fits show that the measured 
peaks are leptocurtic (i.e. are sharper and have longer tails than Gaussian 
distributions), but become very close to a Gaussian distribution for larger 
block size. This confirms the expectation that the macroscopic capillary 
wave model is only valid on length scales larger than some cutoff Bcw 
(compare the discussion after Eq. 0). For a quantitative analysis, a x^-test 
has been performed for the distributions P{h) at temperatures t = —0.05 
and t — —0.01, showing consistency with a Gaussian distribution for block 
sizes B larger than 

Sew ^(5 -7)^, (18) 

which is a reasonable estimate of the capillary wave cutoff Bcw Even for 
large block sizes, however, the value of P{h = 0) is extraordinarily large. 
This is due to the lattice discretization with too low resolution at the center 
of the Gaussian peak. 

The variances of the measured distributions P{h) increase with the 
logarithm of the system size and decrease with the logarithm of the block 
size according to Eq. IZj). The prefactor l/{2TTa) is of the right order 
of magnitude, corresponding to surface tensions a in the range between 
0.033 and 0.038 for t = -0.05. This agrees with the value a = 0.035(1) 
obtained from the scaling formula a = (Jo\t\'^'^ with ag = 1.55(5) [20] and 
V ^ 0.6304(13) E?], and with the value a = 0.0342(2) from ref. 

The capillary wave model thus yields an appropriate description of the 
interface in the three-dimensional Ising model on length scales larger than 
a couple of correlation lengths. 
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4.2 Test of the convolution approximation 

To probe the convolution approximation, Eq. (|5J, the measured interface 
profiles m{z) on various block sizes B and their widths w'^ are analyzed. 
The upper part of Fig. El shows the total interface profiles m{z) (i.e. for 
B = L) for various system sizes L. The profile for L — 32 displays a 
lattice artefact near the interface position, but the profiles for larger sys- 
tem lengths are smoothly varying profiles of a hyperbolic tangent or error 
function type as proposed by the theoretical interface models of Sectional 
The broadening of the profiles with increasing system size due to capillary 
waves is clearly visible. 

The lower part of Fig. |21 shows the interface profiles m(z) on different 
block scales B. In the microscopic regime, i.e. for very small block sizes, 
lattice artefacts dominate and no reasonable interface profile is obtained. 
With increasing block size, however, the profiles approach a smooth hyper- 
bolic tangent or error function type profile which is then broadened with 
increasing block size. This illustrates the renormalization group block spin 
procedure leading from a microscopic Ising "profile" to the mean field pro- 
file, which is then broadened by capillary waves, thus qualitatively following 
the convolution approximation Q. 

For a quantitative analysis the total interface widths li'^ of the total 
interface profiles (i.e. for B = L) are plotted as functions of the logarithm 
of the system size L for various temperatures t, see Fig.|31 As predicted by 
capillary wave theory and scaling, Eq. (|12|l . the width expressed in units 
of twice the correlation length all lie on the same temperature-independent 
straight line, at least for large length scales L. The linear fit displayed in 
Fig.El 

= 0.28(1) log2i + 0.30(4), (19) 

is parallel to the theoretical prediction w'^^, Eq. (|12|l . from the convolution 
approximation with the slope A^w = 0.265(2). 

The logarithmic divergence of the total interface width (for profiles 
with i? = L) in the Ising model has already been observed in other Ising 
Monte Carlo simulations ^1 El , where the temperature closest to the 
critical point was t — —0.2. Converting to temperature-independent units, 
the obtained slopes are icw = 0.222 [El, 0.278 [HI and 0.244 [El. The 
simulation by Stauffer [23 at t = -0.01 and L = 2000 gave the single 
value w'^ — 0.99, which is far below our value = 2.49 (extrapolated 
from Eq. (|19|) '). This discrepancy is due to Stauffer's definition of the 
interface width which only measures the capillary wave part of the width 
and eliminates the intrinsic contribution. 

Fig. ^ shows the widths w'^ of the profiles coarse grained on length 
scales B < L sls a function of the logarithm of the block size 13. For each 
block size B the system shows a width u)^ characteristic for this length 
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scale, independent of the temperature. The variables expressed in units of 
twice the correlation length thus indeed display temperature-independent 
behaviour. Note that the widths of profiles on scales B < L have to be 
treated separately from those of total profiles for B = L due to different 
boundary conditions. For B — L, periodic boundary conditions apply, 
while for B < L the boundary conditions are dictated by the actual system 
configurations, resulting in a different spectrum of capillary waves. 

As expected from capillary wave theory, the widths grow with increasing 
block size B, and for large block sizes a linear increase with log2 B is 
obtained. The linear fit in Fig. 01 

= 0.26(3) log2 5 + 0.57(1) (20) 

is parallel to the theoretical prediction w'^^, Eq. (|12|l . from the convolution 
approximation with the slope A^w = 0.265(2), thus again confirming the 
capillary wave picture on large length scales. 

In comparing the axis intercepts of the linear fits in Figs. El E] to those 
of the theoretical prediction (using as value for the intrinsic width wf^^. 

2 

either the theoretical value j:^' Eq. H12|l . or the "microscopic" value at 
small length scales as determined from Fig. 01 the parameter i?cw can be 
determined. One obtains the reasonable values 

Bcw~(3-7)e (21) 

of the same order of magnitude as in Eq. (|18|l . The cutoff B^w can thus be 
consistently chosen to make theory and simulation agree on large length 
scales. 

4.3 Intrinsic width 

On smaller length scales the width in Fig. ^ grows slower than linear with 
the block size and eventually becomes saturated at a microscopic width 
w 0.2 for microscopic length scales B < This behaviour is in qual- 
itative agreement with the convolution approximation. However, the mi- 
croscopic width is significantly smaller than the intrinsic width, which is 
expected to have the much higher value wf^^ = 7r^/12, see Eq. ((T^ . The 
notion of an intrinsic profile implies that the width should approach this 
value on the intrinsic length scale, i.e. in some region around _B ~ ^. In our 
simulations such an intrinsic regime is not observed, contrary to expecta- 
tions. 

The reason for the deviation from the expected picture is not clear. The 
smallness of the interface width on the smallest scales suggests that it is 
determined by Ising lattice artefacts. This could be clarified by searching 
for the intrinsic profile in the Ising model with larger correlation lengths, 
i.e. closer to the critical point, or in (/'"'^-theory. 
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5 Conclusion 



In this article the interface in the three-dimensional Ising model has been 
studied on different length scales. A procedure has been implemented that 
allows to separate large and small wave length contributions to the inter- 
face profile in a translation invariant way. It is found that on large length 
scales the interfacial properties are well described by capillary wave theory 
with a reasonable choice i3cw ^ (3 — 7) ^ for the cutoff of the capillary wave 
spectrum. The interface width shows universal behaviour, which is in quan- 
titative agreement with the theoretical predictions. The data for the widths 
are consistent with the convolution approximation which includes both the 
intrinsic structure and the capillary waves. The associated intrinsic width 
of the interface, however, turns out to be much smaller than expected from 
field theory, and seems to represent a "microscopic profile" related to the 
discrete nature of the Ising variables. An intrinsic profile could thus not be 
isolated for the Ising model at the simulated temperatures. 

Note added in proof: In ref. ^3 the width of colour flux tubes in the 
3-dimensional Z2 gauge model has been investigated, which is dual to the 
3-dimensional Ising model. The results are in agreement with ours. 
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Figure 1: Distributions P{h) of the interface positions h, coarse grained 
on different block sizes B, in a system of length L = 128 at reduced tem- 
petature t = — 0.01. In addition to the Monte Carlo data (points), Gaussian 
fits are shown (lines). 
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Figure 2: Top: Total interface profiles (i.e. for B = L) for system lengths 
L = 32,64, 128,256 and 512. Bottom: Interface profiles m{z) for different 
coarse graining lengths B = 1, 2, 4, . . . , 128 = L. In both pictures the 
reduced temperature is t = —0.01. 
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Figure 3: Total interface widths ilr^ (in temperature-independent units) 
of profiles for B — L ed, the given temperatures t. Statistical errors are 
smaller than the displayed symbols. For comparison, the intrinsic width 
wf^^ = 7r^/12 and the prediction w'^^^ from the convolution approximation 
(with an offset Acw log2 Bcw) are shown. 
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Figure 4: Interface widths (in temperature-independent units) of pro- 
files coarse grained on block lengths B < L. Symbols are as in Fig.O In 
addition, the correlation length ^ is shown. 
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